Pedestrian index theorem a la Aharonov-Casher 
for bulk threshold modes in corrugated multilayer graphene 
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Zero-modes, their topological degeneracy and relation to index theorems have attracted attention in 
the study of single- and bilayer graphene. For negligible scalar potentials, index theorems can explain 
why the degeneracy of the zero-energy Landau level of a Dirac hamiltonian is not lifted by gauge 
field disorder, for example due to ripples, whereas other Landau levels become broadened by the 
inhomogenous effective magnetic field. That also the bilayer hamiltonian supports such protected 
bulk zero-modes was proved formally by Katsnelson and Prokhorova to hold on a compact manifold 
by using the Atiyah-Singer index theorem. Here we complement and generalize this result in a 
pedestrian way by pointing out that the simple argument by Aharonov and Casher for degenerate 
zero-modes of a Dirac hamiltonian in the infinite plane extends naturally to the multilayer case. 
The degeneracy remains, though at nonzero energy, also in the presence of a gap. These threshold 
modes make the spectrum asymmetric. The rest of the spectrum, however, remains symmetric even 
in arbitrary gauge fields, a fact related to supersymmetry. Possible benefits of this connection are 
discussed. 



I. INTRODUCTION 

Since the experimental realization of graphene it has be- 
come clear that a suspended sheet of graphene is not flat 
but corrugates into a rippled structure. [JQ In the tight 
binding-model for graphene these ripples with their in- 
trinsic curvature lead to a local modification of the hop- 
ping amplitudes. In the low-energy limit given by a Dirac 
hamiltonian the ripples enter as an effective disorder po- 
tential, of which the vector part can be interpreted as 
a nonuniform effective magnetic field. 0] The impact of 
this disorder potential on the spectrum and on transport 
properties has attracted a lot of interest (see ref. 0). For 
example, ref. studies numerically the low-energy spec- 
trum in the presence of certain ripple configurations and 
finds it to be considerably changed when the effective 
magnetic length is comparable to the ripple size. Zero- 
energy Landau-levcl-likc states can then exist within one 
ripple and their degeneracy is no lifted by the inhomo- 
geneity of the effective magnetic field. This should be 
observable as a peak at zero energy in the density of 
states. H, 0] In presence of a scalar potential the degen- 
eracy is lifted, which can also be seen in ref. |^. 

In the quantum Hall problem the small effective mag- 
netic field due to ripples is combined with the strong 
uniform external magnetic field into a total nonuniform 
magnetic field. It is observed that the zero-energy Lan- 
dau level in graphene remains strongly peaked, whereas 
the other Landau levels are broadened, possibly due to 
the inhomogenous field caused by the ripples. 

The stability of zero-energy states — zero-modes — of a 
Dirac electron in a magnetic field of arbitrary shape is 
understood as a consequence of index theorems, which 
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FIG. 1: An imagined typical spectrum for electrons in 
graphene with a gap and a random magnetic field due to rip- 
ples. (The horizontal position carries no meaning.) The two 
valleys K and K' are time-reversed copies in absence of a real 
magnetic field. For ripples without any spatial symmetries 
the spectrum of each valley will in general be non-degenerate 
and of random spacing, except for the threshold modes at 
\E\ = m . Their degeneracy depends only on the total effec- 
tive fiux <1>^ = = 4.1(^0- The rest of the spectrum is 
symmetric around E = due to supersymmetry. 



relate analytical properties of operators to topological 
properties of the space and the fields involved. (See e.g. 
refs. H, [1, [Tol ) There are several index theorems, with 
different range of applicability. Maybe the most famous 
one is the Atiyah-Singer index theorem[T]|. that applies 
to elliptic differential operators on a compact manifold 
of even dimension, e.g. the sphere or the torus. It states 
that for elliptic differential operators with the Fredholm 
property {e.g. Dirac operators Il± = ± iHy on a 
torus, with 11 = — iV + A), the analytical index (the 
number of zero-modes of 11+ , i.e. number of solutions 
to Il^u(x) = 0, minus the number of zero-modes of n_) 
equals the topological index {i.e. the total magnetic flux 
of the gauge field, a topological integer according to the 
Dirac monopole quantization condition on compact man- 
ifolds). In some cases it has been possible to find index 
theorems on non-compact manifolds or for odd space di- 
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mensions. One of many reason for generalizations to 
be interesting is that many Dirac operators occurring in 
quantum mechanics (in particular Dirac operators on the 
infinite plane to be studied here) usually do not have the 
Frcdholm property, aprerequisite for the Atiyah-Singer 
theorem. See refs. 1^, [lO for a review and references on 
these generalizations. 

Fortunately for physicists the essence of this beautiful 
and powerful but rather advanced mathematics is man- 
ifested in some simple examples requiring only elemen- 
tary quantum mechanics, many of them living on infinite 
manifolds. A nice feature of these examples is that the 
wave functions of the zero-modes are obtained explicitly. 
Jackiw and Rebbi[l2| found that in a Id Dirac theory 
there could be a topologically protected zero-mode local- 
ized at a mass soliton. Aharonov and Casherfl^ found a 
very short and simple argument for zero-modes of mass- 
less Dirac electrons in a 2d plane in a nonsingular per- 
pendicular magnetic field of arbitrary shape but finite 
range. The degeneracy of the zero-modes is only deter- 
mined by the total fiux, just like in the Atiyah-Singer 
theorem. In the presence of a mass m these zero-modes 
turn into degenerate threshold modes. Depending on the 
sign of the total flux, they sit either at the E = +m or 
at the E = —m threshold of the gapped spectrum. 

In this paper we focus on Aharonov's and Casher's 
pedestrian argument. In section IIIII we point out how 
simply it also extends to some of the hamiltonians that 
have been considered for multilayer graphene. Stability 
of zero-modes in rippled bilayer graphene was already 
considered by Katsnelson and Prokhorovap^ . there in 
the general but abstract language of the Atiyah-Singer 
theorem, thus applying to a compact manifold, in this 
case the torus resulting from periodic boundary con- 
ditions. Our result is complementary by applying to 
the infinite plane. It also offers a simple generalization 
to multilayers. The pedestrian argument goes beyond 
the Atiyah-Singer theorem by showing that the zero- 
modes are present also when the effective flux though the 
graphene sheet is not an integer. Finally, our argument 
is of pedagogic value as it does not require any higher 
mathematics and also gives concrete wave- functions. 

In section [n] we also extend the discussion by including 
a mass term. In a single layer this corresponds to break- 
ing the sublattice symmetry of the bipartite honeycomb 
lattice, like in experiments on hydrogenated single- 
layers. In a multilayer a gap can also been introduced 
by a transverse potential, like in experiments [l^, on 
gated bilayers. A mass term turns the zero-modes into 
degenerate threshold modes sitting at the gap energy 
(Figs.[T]and[2|). In the quantum Hall problem such a term 
would split the sharp zero-energy Landau-level peak into 
two sharp peaks, symmetrically shifted around zero and 
related to the two valleys (Fig. Because of the val- 
ley degeneracy breaking combination of ripples and an 
external magnetic field the two peaks would be of dif- 
ferent sizes. Such a splitting has actually been observed 
experimentally 0, but has been attributed to other 



mechanisms. 

A known but often not mentioned point is that all the 
mentioned topological arguments for a sharp peak at zero 
energy or at a threshold energy seem to fail if scalar po- 
tentials, for example induced by ripples or by impurities, 
arc not negligible. This important caveat is brought up 
in section [ni but will remain an open question both for 
single-layer and multilayer graphene. 

The third part of the paper makes a note on the sym- 
metry of the spectrum, as illustrated in Figs. [1] and [H 
Apart from the threshold modes, the spectrum of each 
valley remains symmetric around the zero of energy also 
in presence of both a nonuniform vector potential and 
a mass term, provided the mass is constant and that 
the scalar potential is zero. This symmetry in not due 
to the chiral symmetry (7^Hm,Ao,AO-z = -ff_m,-Ao,A 
within each valley, which is only enough to explain this 
symmetry in the massless case for example studied in 
ref. [^. The more general reason of this symmetry we re- 
late instead to the supersymmetric quantum mechanics 
formulated by Wittenp^. Further benefits of this re- 
mark could come from connecting to the rich literature 
on analytic results based on supersymmetry, for exam- 
ple on scattering of Dirac electrons in slowly decreasing 
magnetic fields in which asymptotic states are difficult to 
define. Since the multilayer hamiltonians also have the 
supersymmetric structure, we expect that many such an- 
alytic results should have analogs in the multilayer case. 
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FIG. 2: Same sketch as in Fig. [T] but now including a real 
magnetic field with $em = — 3.2(^o. (The horizontal posi- 
tion only used to guide the eye.) The quasi Landau levels 
are broadened by the ripples, except the Landau level of the 
threshold modes. However, with t&rip — l.l(^o and therefore 

"1>^ = "I-EM -f^rip = -2.1(^0 and = $em - 'I'rip = -4.3(/)o 
the threefold degeneracy of the threshold modes at E = —m 
(E — +m) of valley K (K') becomes twofold (fourfold). The 
rest of the spectrum is in each valley symmetric about E — 
due to supersymmetry. 



II. THRESHOLD MODES IN RIPPLED 
GRAPHENE 

Assume a general magnetic field B{x, y) = d^Ay — dyA^ 
of compact support, i.e. there is a finite disk outside 
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which B{x) = 0. Define the Dirac operators 11 = — zV + 
A and Il± = H^; ± iHy We set ft = e = 1. The argument 
by Aharonov and Casher, that we review in the appendix, 
shows that a massless 2d Dirac hamiltonian 



H = v 



n_ 
n, 



(1) 



has n zero-modes, where n > is only determined by the 
total flux <& = / d^a; B = zt(f>o{n + e) (with = h/e and 
< e < 1) and hence independent of the shape of the 
magnetic field. In the presence of a mass term SH = maz 
these n zero- modes turn into n threshold modes, either 
with energy E = +m or with E = — m, depending of the 
sign of the flux. These solutions are of the form 



u{x) 



u 

^ 




or 



v(x) 





V 





W{x) 



(2) 



(z = {xi -f zx2)/2). / has to be a polynomial for single- 
valuedness and regularity. The magnetic field enters into 
W{x) = J d^x' B{x') In \x — x'\. which far away from 
the region with the magnetic field behaves asymptotically 
as e*^'-^-' ~ la;!^*/"^". Normalizabilty requires the ex- 
ponent to be negative and requires / to be maximally of 
degree n—1. This gives n linearly independent polynomi- 
als and hence n threshold modes, which are zero-modes 
in the massless limit. 

Note that in the case of I'&l/'^iio ~ n + \ — n integer 
flux quanta, the ffi^ solution u oc |a;|^" is not strictly 
square integrable in the plane since the integral diverges 
logarithmically. On a compact manifold this "marginal" 
mode becomes normalizable, therefore the n = <I>/^o-fold 
degeneracy of the lowest Landau level on a torus. 

The valley K of graphene is described by the hamil- 
tonian ([T]) with V Ki lO^m/s, complemented with the 
scalar potential 6H = —IAq and here possibly with 
the addition of a sublattice symmetry breaking term 
5H = maz- The potential m{x) acts with opposite signs 
on the two inequivalent orbitals A and B of the bipar- 
tite honeycomb lattice. The components of the spinor 
(w, v)^ = (V'a' '^b)'^ refer to these orbitals. The gauge 
field {jjL = 0,1,2) is here a sum of the electromag- 
netic part A^,EM and an effective field rip due to 
ripples, [s^l The corresponding hamiltonian for valley K' 
is = —vcr ■ (— iV + Aem — ^rip) — l^o — ™crz acting 
on (^/jg , i/;^ )^ . Thus, the effective magnetic field due 
to ripples changes sign, but not the external magnetic 
field. In absence of an external magnetic field the two 
valleys are mapped onto each other under time-reversal 
ax{H^)*ax = , thus preserving the time- reversal 
symmetry of the total graphene hamiltonian also in pres- 
ence of ripples. With both ^em and Aj-jp nonzero the 
degeneracy between the two valleys is lifted. 

In this context we note that in the absence of a scalar 
potential but with a nonzero gap term maz the thresh- 
old modes of the two valleys sit at the same energy if 



|$rip| > I^emI- For example, with ^lip > the thresh- 
old modes of both valleys sit at £' = —to. In the time- 
reversal invariant case Aem = (Fig. [T|) both valleys 
count the same number of threshold modes and they are 
related by time- reversal V'-,™ = crx{ip-m)*- the op- 
posite case I 'flip I < I'&emI, including the quantum Hall 
scenario, the threshold modes of the two valleys sit at op- 
posite energies (Fig. [5]). Only in the special case <I>iip = 
are they equally many. In general, in a quantum Hall 
problem with a sublattice symmetry breaking gap maz, 
the zero energy Landau-level peak in the density of states 
is split into two symmetrically shifted peaks centered at 
E = +m and E = —m, respectively, but of different 
sizes due to the different fluxes = ^em + ^lip and 
= <i>EM - <i>rip- Experiments [tI. [Tsj flnd at low tem- 
peratures the zero-energy Landau level peak to be split 
into two peaks, resulting in a plateau in the magnetoresis- 
tance at zero doping. We find it suggestive that the two 
peaks appear to be of different sizes. However, the split- 
ting was attributed to counterpropagating edge channels 
dominating the longitudinal resistivitv|2l| or to Zeeman 
splitting in the real spin(20|. 



The effective vector potential derived from the rippling 
is given by Ayip = g2{uxx ~Uyy, 2uxy) [2^ in terms of the 
strain tensor Uij = ^{diUj + djUi + {dih){djh)) with M(a;) 
and h(x) being the in-plane and out-of-plane distortions, 
respectively, g is a coupling parameter that depends on 
the properties of the bonds. Observe that A does not 
change under h —h. A buckled- up ripple gives there- 
fore a flux of the same sign as the equivalent buckled- 
down ripple. Thus, the fluxes of the ripples always add 
up. The flux of one ripple of length I and height h can be 
estimated to be If&l/^o ^ h? /al, with a being the bond 
length of the honeycomb lattice. [2§| Index theorems, or 
the explicit solution found by Aharonov and Casher, ex- 
plain why ripples can lead to degenerate zero-modes or 
threshold modes. 

The accompanying effective scalar potential 5H = 



igi{ux 



y) destroys this exact degeneracy (as it is 



nonuniform) . An approximate degeneracy could still hold 
if gi is small enough. However, it is far from obvious 
that this is the case for monolayer graphene, considering 
estimates (2^ for carbon nanotubes with gi « 30 eV and 
g2 ~ 1.5 eV. On the other hand, it is not impossible that 
the bare value of the scalar potential could be substan- 
tially reduced for example through screening. As far as 
we know this has not been clarifled. It remains therefore 
unclear to us if the observed sharpness of the zero energy 
Landau level is really due to index theorems, with some 
mechanism suppressing the scalar potential, or if it must 
be attributed to some other reason. As for multilayers, 
to which we now turn to, we do not know of any such 
estimates. With the above caveat mentioned, we now go 
on to analyze threshold modes in multilayers for the case 
that the scalar potentials would turn out to be negligible 
there. 
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III. THRESHOLD MODES IN MULTILAYER 
GRAPHENE 

The intense study of monolayer graphene has been fol- 
lowed by that of tight-binding models for bilayer and 
multilayer graphene, with low-e nergy continuum models 
similar to the monolayer one. [2^ [25n26l. \2l\ These multi- 
layer sheets have exotic properties of their own, but offer 
also an interesting interpolation back to the source ma- 
terial graphite. In the simplest tight-binding model, the 
electrons at the valley K in a bilayer are in a low-energy 
range described by the effective hamiltonian 

with J = 2.(2^ This hamiltonian acts on the spinor 
(V'AjV'b) i'^ ^ bilayer arranged so that the B-orbitals of 
the AB-layer are placed on top of the A-orbitals in the 
AB-layer, giving the leading inter-layer tunneling chan- 
nel. The superposed orbitals A and B dimerize in the 
low-energy limit, leaving effectively the A and B orbitals. 

For a three layers and beyond [H, [H, H^l there are sev- 
eral ways of stacking the layers, for example the Bernal 
stacking (ababa...) — the most common form in natural 
graphite — or the rhombohedral stacking (abcabc...), with 
a, b, and c denoting the three different placements of 
honeycomb sheets that can occur in a stacking. With 
the simplest approximation for inter-layer tunneling one 
obtains for a rhombohedrally stacked A^-layer the effec- 
tive hamiltonian ([3]) with J = N. For a Bernal stacked 
TV-layer one finds instead 

H,s - J, , = ^ (4) 

Ji i 

with Ji — 2, except Ji = 1 if iV is odd — thus a tensor 
product of bilayer and monolayer hamiltonians. Other 
stackings give a structure that is intermediate between 
the rhombohedral and the Bernal one.p7| 

Also multilayer graphene sheets form ripples [2^, and 
the existence and stability of zero-modes is interest ing 
to address. Particularly, Katsnelson and Prokhorova|14j 
gave a formal proof, based on the Atiyah-Singer index 
theorem complemented with a result [29| on indices of 
composed elliptical operators. They showed that the 
mentioned bilayer hamiltonian on a compact 2d manifold 
(in particular the torus obtained from assuming periodic 
boundary conditions) has zero-modes, precisely twice as 
many as the monolayer hamiltonian ([1]). The contribu- 
tion of the present paper is to note in a pedestrian way 
(although now for an infinite plane) that this fact and the 
generalization to arbitrary J follow from the argument 
of Aharonov and Casher with the following straightfor- 
ward extension. Assume <& to be negative. There are 
the n zero-modes (u, v)^ = {f{z)e^^'^\0)'^ that satisfy 
n+/e^ = 0. Since n+ = -id;,. + A+ it follows that 

n+{z*yfe'^ = -ij(zy~^fe'^. (5) 



As a consequence, all j = 0, . . . , J— 1 give wave functions 
u = (z*)-'/(z)e^(^) that are zero-modes of Il-^u = 0. 
The hamiltonian ^ has therefore nJ zero-modes and 
the total hamiltonian ^ has nN zero-modes, indepen- 
dently of the shape of the gauge field and independent 
of the stacking configuration. In the presence of a con- 
stant mass term maz in the hamiltonian, these zero- 
modes turn into degenerate threshold modes, exactly as 
for monolayer graphene. 

In the special case of a constant [13 negative magnetic 
field B the subspace of zero-modes is nothing but the 
J first Landau levels with 11+ /\/—2B being a lowering 
operator of the Landau level index. The Dirac Landau 
levels with energy E oc iy^ 7^ are given by (u, u)"^ = 
{(pp zt(pj_j)'^ with j > J and in terms of the Landau 
level j wave-functions (t)j{z, z*) of a spinless Schrodinger 
hamiltonian. The zero-modes are given by (v^j, 0)'^ with 
j = 0, . . . , J — 1. In particular, the zero energy LL of 
a bilayer has twice the degeneracy of the corresponding 
monolayer analog. [2^ 

IV. A SUPERSYMMETRIC SPECTRUM 

In addition to the note on threshold modes in multilayers, 
we will also make a note on the rest of the spectrum. The 
chiral symmetry 

(^zHm,Ao,A(^z = —H^„i,~Aq,A (6) 

within each valley and valid for all J implies that the 
symmetry of the spectrum around i? = is broken by a 
scalar potential or a mass term. (For the massive case 
and J odd, the particle-hole conjugation a^H^ A^'^^ ~ 
—Hm.-Ai^ guarantees a symmetric spectrum per valley, 
if instead gauge potentials are absent. For J even, let 
<7x "* <^y) In particular, there are threshold modes either 
only at E = -|-to or only at _B = — m, depending on the 
sign of the total flux. 

Amazingly, however, the rest of the spectrum within 
each valley remains symmetric, even when m 7^ 0, pro- 
vided the mass term m is constant and the scalar poten- 
tial ^0 is zero. With these conditions the property 

0=[H\az] = {H,^[H,az]}, (7) 

is fulfilled, which we note holds for any J and arbitrary 
A. It implies that the hermitian operator 

-n-) (8) 

maps positive energy states into negative energy states, 
except the threshold modes, which it kills. In the mass- 
less limit it coincides with the chiral symmetry ([S]), ex- 
cept for zero-modes. 

To be more fancy, the symmetry is a manifestation of 
supersymmetry as defined in supersymmetric quantum 
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mechani cs. [Tgll The Dirac hamiltonian is one important 
example[30l, but we can obviously generalize to the mul- 
tilayer hamiltonians. This has already been noted [sH in 
the study of the quantum Hall spectrum in multilayers. 
However, the quantum Hall spectrum is highly degener- 
ate, and we want to stress that the supersymmetry in 
multilayers holds for arbitrary magnetic fields and that 
the minimum message of supersymmetry is most clearly 
seen in random fields. The real usefulness of redressing 
([7]) into supersymmetry would be if one could make con- 
tact with the rich literature on analytic results based on 
the latter. (Supersymmetry, for instance, explains why 
the Dirac equation for some potentials, in particular 3d 
Dirac electrons in a Coulomb potential, can be solved ex- 
actly and the spectrum can be constructed algebraically.) 
One thing that might be interesting for the study of rip- 
ples could be results on scattering of Dirac particles in 
slowly decreasing magnetic fields when asymptotic states 
are not easy to define, see e.g. ref.[l3. For results relying 
on supersymmetry, we expect that similar results should 
hold for the supcrsymmctric hamiltonians of multilayer 
graphenc. 

The essential structure is that there is a unitary sclf- 
adjoint operator r (in our case a^) with = 1 for 
which there is a decomposition of H in hermitian parts 

H = iJodd + ^^ovon such that [i/cvon,^] = {-ffodd, t} = 

{i?evon, ^^odd} = 0. In such a case one can go on to form 
the supersymmetric hamiltonian H. = ^H^^^ = = 
2Qi = {2,2'''} and the supercharges 2i = ^-^odd, 
22 = i[Qi,T] and 2 = i(2i +«22) with {2i,22} = 
and ~ 0. Wc recognize 22i as the massless Dirac 
hamiltonian and 22 as the operator ([5|) conjugating the 
spectrum of the massive Dirac hamiltonian. 2 and 2^ 
are fermionic ladder operators to be discussed below. For 
the hamiltonian H = Dx{x)ax+Dy{x)ay+maz (for mul- 
tilayers D± = one finds 







2T-L 



(9) 



A supersymmetric hamiltonian 7Y is obviously positive 
definite. Because of 2] = 0, all positive energy eigen- 
values Ei of Ti correspond to a degenerate doublet 



(10) 



one state coined "bosonic" and the other "fermionic". 
The supercharges 



2 = 









and 2 



t _ 







(11) 



act as fermionic ladder operators stepping between them. 
The two states in one doublet are linear combina- 
tions of two eigenstates of H related by 22- (In par- 
ticular, in the case of constant B the doublets are 
{(^^, 0)T, (0, 

The exceptions are the zero-modes of Ti, also zero- 
modes of all the supercharges. (Thus, the zero-modes 



of Hm=o or the threshold modes of H.) Supersymmetry 
does not imply their existence. If they exist they do not 
need to come in pairs since they are killed both by 2 and 
by 2^. (Equivalently, the threshold modes of H cannot 
be conjugated by 22-) 

The asymmetry in number of bosonic and fermionic 
zero-modes is the Witten index, which equals the Atiyah- 
Singer index for the operators D+ when they have the 
Fredholm property. {Sec e.g. ref.uQ.) For Dirac opera- 
tors on an infinite plane, the Witten index defined as the 
asymmetry between bosonic and fermionic states should 
equal the number of zero-modes given by the Aharonov 
and Casher argument. However, with the analytical defi- 
nition of the Witten index discussed in ref. [13, it actually 
remains ^/4>q also on the infinite plane. The Witten in- 
dex is then only related but not identical to the number 
of zero-modes according to Aharonov and Casher. 



V. SUMMARY 

In this paper we pointed out that Aharonov-Casher argu- 
ment for zero-modes for 2d Dirac electrons in a magnetic 
field generalizes naturally to the low-energy hamiltonians 
studied in the case of multilayer graphcne. We also dis- 
cussed the relationship of this work to that of Katsnelson 
and Prokhorovafl^ . who gave a formal proof for bilaycrs 
based on the Atiyah-Singer index theorem. We extended 
the discussion to the presence of a uniform gap, in which 
case the degeneracy remains and might still be observable 
as a peak in the density of states, but at nonzero energy 
with the zero- modes instead being threshold modes. 

Further, we made a note on the symmetry of the spec- 
trum within each valley in the presence of an arbitrary 
magnetic field. In the massless case the chiral symmetry 
guarantees a symmetric spectrum. Apart from the asym- 
metry of the threshold modes, the multilayer spectrum 
remains symmetric for a uniform nonzero mass term. We 
related this to the supersymmetric structure of the con- 
sidered multilayer hamiltonians. We expect that many 
of the analytic results based on the supersymmetry of 
the Dirac hamiltonian should have analogs applying to 
the hamiltonians of multilayer graphcne. Interesting ap- 
plications might be found for example for the study of 
scattering in a nonuniform magnetic background. 

All these arguments fail if the scalar potential is not 
negligible compared to the vector potential. We brought 
up this issue but do not know of any solid answers neither 
for single-layer nor for multilayer graphenc. 
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author is grateful to E. Mariani for discussions about the 
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Appendix: Aharonov's and Casher's argument 

The two-dimensional Dirac equation = (H — E)ip for 
electrons of mass m and charge —1 is determined by the 
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hamiltonian H = UxO'x + ^yO'y + maz- {c = h = e = 1, 
hence 4>o = h/e = 2tt.) The threshold modes at E = +m 
have to satisfy ip+m ~ {u, 0)^ and H+u = 0. Like- 
wise, there are threshold modes E = — m of the form 
tp-m = {0,v)'^ provided n_w = 0. It will now be in- 
vestigated which of these possibilities really gives nor- 
malizable wave fmrctions. Introduce z = {x + iy)/2, i.e. 

= dzdz' and 11+ = {-id^* +A+) with A± = Ax±iAy. 
The Ansatz u{x) = f{z)e^^^\ with f{z) an arbitrary an- 
alytic function, satisfies = II+m = fe^{—idz*W + A+). 
Acting with on d^* W = —iAj^ gives the Poisson equa- 
tion V^W ^ B — iV ■ A. Thanks to the boundary condi- 
tions given by the compact support of B this equation can 
be inverted with the help of In |a; — a;'| = 2tt6{x — x'), 
resulting in 

W{x) =F{z) + Giz*) + 

+ {B{x')-iW ■A{x'))ln\x-x'\, ^ ' 

with F[z) and G{z*) arbitrary analytic functions. We 



put F{z) = as it is already accounted for by f{z). 
Also, dz'W ~ —iA+ implies G{z*) = 0. Similarly, the 
Ansatz v{x) = f{z*)e~^ ^^'> results in the same W as 
for u. Choosing Coulomb gauge V • A = implies W* = 
ly = ^ / (Px'B{x') In \x - x'\. 

Far away from region of flux (i.e. \x\ > \x'\ and 
ln|a; — x'\ ~ hi|a;|) one has asymptotically e^^*^^^ ~ 
|^|±*/2ir^ Normalizability imposes that only E = +m 
[E = —to) threshold modes can exist for $ < (<& > 0). 
The degeneracy of these modes comes from the possi- 
ble choices of f{z). Single-valucdness of the wave func- 
tion requires f{z) = X^sez'^s^^" and normalizability at 
oo requires degf < \^\/(j)o — 1 = n + e — 1, where 
< e < 1. Normalizability within the region of flux re- 
quires B(x) to be non-singular and all powers in / to 
be non- negative. Thus, / is a polynomial in z (or z*) of 
maximal degree ?i — 1. There are n linearly independent 
polynomial dcg/ < ^ ~ 1, therefore the n-dimensional 
subspace of threshold modes. 
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